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We present a new method for calculating the bubble nucleation rate in first order 
phase transitions non-perturbatively on the lattice. The method takes into account 
all fluctuations and the full dynamical pre-factor. We also present results from 
applying it to the cubic anisotropy model, which has a radiatively induced, strongly 
first order phase transition. 



Electroweak baryogenesis in the Standard Model (or extensions of it) happens 
on or near the surface of the bubbles which nucleate and subsequently grow 
during the first order Electroweak phase transition!!! In this case, it is very 
difficult to compute the bubble nucleation rate analytically with sufficient ac- 
curacy: since the transition is radiatively generated, no classical bubble so- 
lution exists, and the whole idea of separating the classical bubble from the 
fluctuation determinant, needed for Langer's nucleation theoryp becomes cum- 
bersome. Furthermore, the long-distance physics of the Electroweak theory is 
inherently non-perturbative, making lattice simulations necessary. 

How can one calculate the nucleation rate on the lattice? The most 
straightforward method is to take an ensemble of configurations in the meta- 
stable state, evolve each with Hamilton's equations of motion, and wait for 
tunneling to happen. This has been done, for instance in the fsing modelta and 
the cj) 4 models However, the cooling rate of the Universe during the Electroweak 
phase transition is many orders of magnitude smaller than the timescale of mi- 
croscopic interactions. Thus, the system has plenty of time for probing its 
phase space, and the tunneling will happen through very strongly suppressed 
configurations (p ~ e -100 ) after a small amount of supercooling. However, it 
is impossible to make Monte Carlo simulations with e 100 iterations! 

Below, we shall describe a Monte Carlo method which fully overcomes this 
problem, and which can be used to calculate both the static and the dynamical 
parts of the nucleation rate. Instead of using the full SU(2) x f/(l)+Higgs, 
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we chose as our toy model the cubic anisotropy model, which has the desired 
features: 

• It has a radiatively induced, first order transition. As mentioned above, 
this makes analytical computations quite difficult, see refta 

• It is simple to simulate on the lattice. 

The continuum Hamiltonian of the cubic anisotropy model in three dimensions 
is {4>\ and <j>2 are scalar fields, ~k\ and 7T2 the associated canonical momenta): 

0=1,2 

The Hamiltonian is discretized on the lattice using an improved (next-to- 
nearest neighbour) Laplacian. Included also are loop corrections to second 
order, which enter in the lattice versions of Ai, A2, <pi, 4>2 and m 2 . We fix 
A2 = 8A1, giving us a strongly first order phase transition between the sym- 
metric and broken phases, as we vary m 2 . In this case m 2 plays the role of the 
temperature. We find the critical m 2 , m 2 , the latent heat_and the interface 
tension. The full results will be presented in a future papem 

The calculation of the rate can be split up into two main parts, the static 
probability of the critical bubbles (part /), and the dynamical flux of the phase 
space through the bubble configurations (part II). As an order parameter, we 
choose the space average of $ 2 = <f> 2 + §\. 

• The non-dynamical part I is the probability of having a critical bubble con- 
figuration. It is found by calculating the probability distribution P($ 2 ) of 
the order parameter (see Fig. |l|), and defining the critical bubble to be at 
the minimum, Z . We define a region [Z — |, Z + |] as the critical bubble 
region. 

At the critical temperature, because of the finiteness of the lattice, the most 
suppressed configurations are slabs (see Fig. [j]), and the bubble configura- 
tions are found on the inner sides of the two peaks. Lowering the temperature 
(m 2 ), the critical bubble becomes smaller and eventually it will fit inside the 
lattice volume. Part I is then the probability density in the critical bubble 
region (=-*■) divided by the probability of the metastable state (Pa)- 

I = , Pz= f +2 F($ 2 )d$ 2 , P A = f P($ 2 )d$ 2 (2) 

PAe Jz-i Jo 

P(4> 2 ) is found by using multicanonical methods! It can be reweighted to 
several different values of m 2 , and part J(m 2 ) calculated. 
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Figure 1: Sketch of the different geometries of the configurations in a finite box (left). The 
probability distribution is reweighted to several different values of m 2 where the critical 
bubble fits inside the lattice volume. For each m 2 , Z is determined and e is chosen (right). 



• Part Ila is the flux of order parameter through the critical bubble value, Z. 
It can be calculated analyticallp: 

IIa=-{\ — \) = J— 3 
2 at V 7tV 

• However, not all of the flux in Ila leads to tunneling. The bubble wobbles 
back and forth for some time before falling into one of the minima. Since 
we are actually sampling the configurations at <I> 2 = Z, we are over-counting 
the true phase space flux from the symmetric phase to the broken phase by 
the average number of crossings of Z the trajectories make. Therefore, we 
use the Hamiltonian in Eq. (]]]) to calculate the real time trajectories forward 
and backwards in time from initial configurations in the critical bubble region 
[Z — |, Z + |]. Gluing the two together into a full trajectory, we can then 
count the number of crossings, and determine whether the trajectory starts 
and ends in different minima. Finally, we can correct the part Ila with part 
lib: 

lib ee iV^ (4) 

where N is the number of trajectories, Ci is the number of crossings of the 
i'th trajectory and Si is 1 if the trajectory tunneled, if not. 
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Real-time trajectories. X 2 a = 3, m L = 0.081 . 
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Figure 2: From the trajectories we can count the crossings and see if it is a tunneling 
trajectory. Without the noise term, the bubble stays on the barrier too long and the number 
of crossings becomes too large. 



However, since the Hamiltonian evolution of the bubble conserves the total 
energy, the finite size of the system poses a problem: when the bubble grows 
(shrinks), it releases (absorbs) latent heat, and the temperature grows (de- 
creases) . This effect tends to unphysically "stabilize" the bubble around the 
critical radius, see Fig. |^. This can be solved by going to very large volumes, 
which is too expensive, or by adding a small amount of thermal noise in the 
system. The thermal noise consists of an update of the momenta: 



tt' = + e£ , e=J=^ (5) 

V -L/ 

where £ is taken from a Gaussian of width 1. The amplitude of the noise 
is tuned so that the tt's are thermalized after evolution of length t = X^aL. 
The final rate is insensitive to the precise value of e. 



The final rate is (see Figure ||): 



_ / x 11a x lib , . 

r = v (6) 
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Figure 3: The full nucleation rate compared to the thin- wall approximation (line). 



In Fig. y we compare our result with the thin- wall approximation. At smallest 
supercoolings we have studied, the thin-wall result gives a rate ~ 25 orders of 
magnitude too large. However, for a fixed rate, the thin-wall approximation 
gives "only" ~ 30% too small supercooling Sm 2 . 

To conclude, we have described a lattice Monte Carlo method for calcu- 
lating the full bubble nucleation rate in first order phase transitions. The 
method is especially well suited for transitions which happen through very 
strongly suppressed bubble configurations; indeed, most first order phase tran- 
sitions in Nature fall into this class. We have applied this method to the cubic 
anisotxopy model in 3D, for a full description we refer to the forthcoming 
paperEl This method has also been applied to SU (2)+Higgs theory, hut with 
a heat bath evolution instead of the Hamiltonian real time evolutionQ 
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